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BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


1  Scientific  Progress  and  Accomplishments 

In  this  program  two  major  experimental  goals  and  one  important  theoretical  study 
have  been  accomplished  in  the  past  three  years.  The  important  experimental  results 
include  the  first  studies  of  direct  evaporative  cooling  of  fermions  in  an  all-optical  trap 
and  the  achievement  of  all-optical  production  of  a  degenerate  Fermi  gas.  In  addition, 
the  first  scaling-law  theory  of  evaporation  in  time-dependent  optical  traps  has  been 
developed. 

These  achievements  pave  the  way  for  important  studies  of  trapped  Fermi  gas 
analogs  of  high  temperature  superconductivity  in  condensed  matter  systems.  Since 
the  density,  temperature,  and  interaction  strength  are  widely  tunable  in  atomic  sys¬ 
tems,  pairing  interactions  spanning  the  range  from  weak  to  resonantly  enhanced  will 
be  possible.  Further,  trapped  multi-component  Fermi  gases  can  be  used  to  model  few 
body  interactions  in  nuclear  matter.  Such  systems  will  provide  widely  variable  data 
for  comparision  with  new  nonperturbative  calculation  methods  recently  developed  for 
few-body  scattering  in  nuclear  and  high  energy  physics. 

1.1  Overview 

This  program  explores  fundamental  interactions  in  multi-component  Fermi  gases. 
Since  single  component  Fermi  gases  do  not  interact  at  low  temperatures  as  a  conse¬ 
quence  of  the  Pauli  exclusion  principle,  it  is  necessary  to  trap  and  cool  more  than 
one  fermionic  species  in  order  to  study  interactions.  In  current  experiments,  40K 
and  6Li  are  considered  particularly  promising  for  observation  of  superfluidity  in  a 
Fermi  gase,  because  these  atoms  exhibit  magnetically  tunable  Feshbach  resonances  in 
their  scattering  interactions.  The  predicted  superfluidity  is  an  atomic  gas  analog  of  a 
superconducting  transition.  For  systems  with  resonant  interactions,  transition  tem¬ 
peratures  of  up  to  half  the  Fermi  temperature  are  predicted.  This  would  make  these 
atomic  systems  the  highest  temperature  superconductors  ever  studied,  in  units  of 
the  Fermi  temperature.  However,  the  most  interesting  mixtures  of  these  Fermi  gases 
must  be  confined  in  an  all-optical  trap,  because  the  states  of  interest  are  repelled 
from  magnetic  traps. 

In  the  previous  program,  we  developed  the  first  stable  all-optical  trap  for  neutral 
atoms,  which  is  based  on  a  stable  CO2  laser  [1],  Using  this  trap,  a  gas  of  fermionic 
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6Li  has  been  confined  with  a  1/e  lifetime  of  400  seconds.  Building  on  this  result,  a 
primary  goal  of  the  current  program  became  the  development  of  evaporative  cooling 
methods  which  could  be  applied  directly  in  the  stable  CO2  laser  trap.  This  effort 
culminated  in  the  first  all-optical  production  of  a  degenerate  Fermi  gas  last  year. 

The  spectacular  result  of  all-optical  production  of  a  degenerate  Fermi  gas  has 
been  highlighted  in  several  places,  including  the  AIP  Physics  News  Update  (March 
13,  2002),  Physics  Today  (May,  2002),  and  in  the  Photonics  News  section  of  Photonics 
Spectra  (May,  2002).  In  addition,  the  work  has  been  highlighted  in  longer  article  in 
Photonics  Spectra  (June,  2002). 

In  the  following,  the  primary  results  of  the  current  program  are  outlined  briefly. 

1.2  Scaling  Laws  for  Evaporative  cooling  in  Time-dependent 
All-Optical  Traps. 

In  the  past,  the  approach  to  preparing  degenerate,  optically  trapped  degenerate  gases 
has  been  to  achieve  degeneracy  in  a  magnetic  trap,  and  then  to  transfer  the  gas  to 
an  optical  trap.  Subsequent  state  preparation  is  accomplished  by  employing  radio¬ 
frequency  transitions  to  produce  the  required  mixture  of  atomic  states.  In  this 
way,  states  which  cannot  be  confined  in  a  magnetic  trap  are  produced  in  a  state- 
independent  all-optical  trap.  At  present,  no  group  has  been  able  to  prepare  a  system 
suitable  for  study  of  Fermi  superfluidity  by  this  method.  A  key  problem  is  heating 
that  arises  during  the  transfer  from  the  magnetic  trap  to  the  optical  trap.  In  ad¬ 
dition,  imperfections  in  state  preparation  can  lead  to  inelastic  decay  of  some  atoms 
and  subsequent  heating.  This  cumbersome  procedure  can  be  greatly  simplified  by 
directly  producing  the  required  degenerate  mixtures  by  evaporation  in  an  optical 
trap.  Further,  the  all-optical  method  enables  trapping  and  cooling  of  systems  such 
as  diamagnetic  molecules  which  cannot  be  magnetically  confined. 

We  developed  the  first  scaling  law  theory  of  evaporation  in  all-optical  traps  which 
are  adiabatically  lowered  during  the  evaporation  sequence.  Unlike  the  scaling  laws 
developed  previously  for  evaporation  in  constant  depth  magnetic  traps,  we  include 
the  loss  of  energy  arising  both  from  evaporation  and  from  energy  transfer  to  the 
trapping  potential.  The  latter  assures  that  the  phase  space  density  remains  constant 
when  the  trap  is  lowered  without  evaporation.  The  predictions  show  that  evaporation 
in  all-optical  traps  can  be  very  efficient,  leading  to  degeneracy  with  relatively  little 
atom  loss. 

For  atomic  systems  with  a  moderately  high  collision  rate  and  a  high  initial  phase 
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space  density,  the  scaling  laws  predict  that  degeneracy  should  be  readily  attainable 
by  lowering  the  trap  depth  by  a  factor  of  about  100.  These  predictions  were  made 
as  part  of  Ken  O’Hara’s  thesis  a  year  prior  to  the  demonstration  of  a  BEC  in  an 
all-optical  trap  by  the  Georgia  Tech  group,  who  used  evaporative  cooling  in  a  trap 
lowered  by  about  a  factor  of  100.  Our  predictions  are  consistent  with  the  results  of 
the  Georgia  Tech  group.  A  Rapid  Communication  describing  the  scaling  law  model 
was  published  last  year  by  our  group  [2]. 

1.3  All-Optical  Production  of  a  Degenerate  Fermi  Gas. 

Last  year,  we  succeeded  in  preparing  a  degenerate  two-component  mixture  of  6Li 
fermions  by  direct  evaporation  in  a  stable  CO2  laser  trap  [3].  We  trap  the  two 
lowest  hyperfine  states  which  are  predicted  to  have  a  magnetically  tunable  Feshbach 
resonance  in  the  scattering  length.  The  trapped  mixture  has  the  property  that  the 
scattering  length  is  zero  at  zero  magnetic  field,  enabling  interactions  to  be  turned 
on  and  off  at  will.  This  system  is  ideally  suited  for  study  of  recently  predicted 
resonance  superfluidity  and  for  study  of  the  transition  from  weak  BCS  pairing  to 
Bose  condensation  of  strongly  bound  pairs. 

In  the  experiments,  we  are  able  to  trap  up  to  3.5  x  106  atoms  in  the  CO2  laser 
trap  by  loading  from  a  6Li  magneto-optical  trap.  Free  evaporation  is  initiated  by 
applying  a  100  G  bias  field.  After  6  seconds  of  evaporation  at  fixed  trap  depth,  we 
obtain  1.3  x  106  atoms  at  a  temperature  of  50  and  an  initial  phase  space  density  of 
8  x  10-3  per  state  which  is  very  high.  Forced  evaporation  is  accomplished  by  lowering 
the  trap  laser  intensity  over  50  seconds.  Using  time-of-flight  imaging,  we  find  the 
temperature  drops  to  15  after  10  seconds  of  evaporation  yielding  T/TF  =  1,  where 
Tf  is  the  Fermi  temperature.  After  40  seconds  of  forced  evaporation,  T/TF  <  0.5 
with  105  atoms  remaining. 

1.4  Scattering  Interactions  at  High  Magnetic  Fields 

Recently,  we  installed  a  new  high-field  magnet  system,  which  produces  a  uniform 
field  of  up  to  1100  G  in  the  trap  region.  Using  this  system,  we  are  exploring  elastic 
scattering  and  loss  in  our  two-state  Fermi  gas.  In  the  field  region  near  650  G,  there  is 
loss  on  a  time  scale  of  a  second  at  temperatures  of  5  /jK  and  densities  of  3  x  1013/cm3. 
However,  at  fields  just  above  the  predicted  Feshbach  resonance  at  850  G,  there  is 
negligible  loss  and  heating.  This  bodes  well  for  experiments  on  superfluidity. 
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We  have  been  able  to  evaporate  at  a  field  of  900  G,  achieving  degeneracy  in  just 
a  few  seconds  as  a  result  of  the  very  large  elastic  scattering  cross  section  at  this  field. 
We  achieve  temperatures  less  than  3  fiK  in  a  50-50  mixture  of  the  two  lowest  hyperfine 
states,  where  the  Fermi  temperature  is  more  than  9  /dK.  These  experiments  appear 
close  to  producing  the  requirements  for  observing  superfluidity  in  a  Fermi  gas. 
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2  Invited  Talks  and  Seminars 

During  this  period,  our  group  presented  9  invited  talks  and  seminars  discussing  this 
research: 

1)  “Quest  for  superfluidity  in  an  optically  trapped  Fermi  gas,”  Seminar  on  Modern 
Optics  and  Spectroscopy,  M.  I.  T.,  (Cambridge,  MA  September  19,  2000). 

2)  J.  E.  Thomas,  “Optically  trapped  Fermi  gas,”  Institute  for  Theoretical  Atomic 
and  Molecular  Physics  and  Physics  Department,  Harvard  University  (Cambridge, 
MA  September  20,  2000). 

3)  J.  E.  Thomas,  “Optical  Trapping  of  a  Two-Component  Fermi  Gas,”  Optical  Physics 
Interdisciplinary  Laboratories  Seminar,  (U.  Michigan,  April  3,  2001). 

4)  J.  E.  Thomas,  S.  R.  Granade,  M.  E.  Gehm,  M.-S.  Chang,  and  K.  M.  O’Hara,  “Op¬ 
tical  Trapping  of  a  Two-Component  Fermi  Gas,”  Fifteenth  International  Conference 
on  Laser  Spectroscopy  (Snowbird,  Utah,  June  10-16,  2001),  invited  talk. 

5)  K.  M.  O’Hara,  M.  E.  Gehm,  S.  R.  Granade,  M.-S.  Chang,  and  J.  E.  Thomas, 
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“Coherence  in  an  optically  trapped  Fermi  Gas,”  Eighth  Rochester  Conference  on  Co¬ 
herence  and  Quantum  Optics  (Rochester,  NY,  June  13-16  2001),  invited  talk. 

6)  K.  M.  O’Hara,  M.  E.  Gehm,  S.  R.  Granade,  and  J.  E.  Thomas,  “Two-component 
Fermi  gas  in  an  optical  trap,”  17th  Interdisciplinary  Laser  Science  Conference,  (Long- 
beach,  CA  October  14-18  2001),  invited  talk. 

7)  J.  E.  Thomas,  “Evaporative  cooling  of  a  two-component  6Li  Fermi  gas  in  a  stable 
CO2  laser  trap,”  Joint  Workshop  hosted  by  the  Institute  for  Theoretical  Atomic  and 
Molecular  Physics  and  the  Harvard-MIT  Center  for  Ultracold  Atoms:  Beyond  BEC: 
Ultracold  atoms  beyond  mean-field  physics,  (Harvard-Smithsonian  Center  for  Astro¬ 
physics,  November  2-3,  2001),  invited  panelist. 

8)  J.  E.  Thomas,  “The  quest  for  superfluidity  in  a  trapped  Fermi  gas,”  (Argonne 
National  Laboratory,  November  5,  2001). 

9)  J.  E.  Thomas,  “The  quest  for  Fermi  superfluidity:  Implications  of  an  all-optical 
degenerate  Fermi  gas,”  DAMOP  2002  (Williamsburg,  May  2002),  invited  talk. 

3  Technological  Applications 

Our  current  experiments  explore  the  dynamics  of  an  ultracold,  trapped  fermionic 
vapor  of  6Li.  Potentially,  these  studies  will  enable  investigation  of  superfluid  pairing 
in  a  two-component  vapor  with  tunable  interactions.  6Li  has  recently  been  shown  to 
be  the  highest  temperature  superconductor  known,  when  the  transition  temperature 
is  measured  in  units  of  the  Fermi  temperature.  Hence,  studies  in  this  system  are 
likely  to  lead  to  new  and  fundamental  insights  into  the  nature  of  high  temperature 
superconductivity  which  strongly  impacts  materials  science.  Further,  the  ultrastable 
CO2  laser  enables  trapping  of  both  atoms  and  molecules  with  broad  applications  to 
precision  measurements  and  novel  clocks.  For  example,  new  clocks  can  be  based 
on  coherent  superposition  states  of  fermions  which  are  prevented  from  colliding  by 
the  exclusion  principle.  Our  experiments  are  also  well  suited  for  exploring  novel 
matter-wave-optical  processes,  such  as  abrupt  transitions  between  interacting  and 
noninteracting  states  on  a  time  scale  short  compared  to  the  time  over  which  an  atom 
travels  a  deBroglie  wavelength.  In  this  case,  novel  many  body  quantum  dynamics 
is  expected.  Such  experiments  will  enable  unique  studies  of  nonlinear  atom-wave 
processes,  enabling  new  techniques  for  manipulation  and  control  of  matter-wave  fields 
with  potential  applications  in  nanolithography. 
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All-Optical  Production  of  a  Degenerate  Fermi  Gas 
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We  achieve  degeneracy  in  a  mixture  of  the  two  lowest  hyperfine  states  of  6Li  by  direct  evaporation  in 
a  CO2  laser  trap,  yielding  the  first  all  optically  produced  degenerate  Fermi  gas.  More  than  105  atoms  are 
confined  at  temperatures  below  4  pK  at  full  trap  depth,  where  the  Fermi  temperature  for  each  state  is 
8  /iK.  This  degenerate  two-component  mixture  is  ideal  for  exploring  mechanisms  of  superconductivity 
ranging  from  Cooper  pairing  to  Bose-Einstein  condensation  of  strongly  bound  pairs. 


DOI:  10. 1 1  Q3/PhysRevLett.88. 120405 

Degenerate  two-component  Fermi  gases  offer  tantaliz¬ 
ing  possibilities  for  precision  studies  of  pairing  interactions 
in  systems  for  which  the  density,  temperature,  and  interac¬ 
tion  strength  are  widely  variable.  Of  particular  interest  are 
certain  two-component  mixtures  of  and  6Li  which  ex¬ 
hibit  magnetically  tunable  Feshbach  resonances,  enabling 
variation  of  the  s-wave  scattering  interaction  from  strongly 
repulsive  to  strongly  attractive.  Attractive  mixtures  in  these 
systems  are  analogs  of  superconductors,  since  they  have 
been  predicted  to  undergo  a  superfluid  transition  as  a  re¬ 
sult  of  Cooper  pairing  at  experimentally  accessible  tem¬ 
peratures  [1,2].  Recently,  two  groups  have  predicted  the 
possibility  of  superfluidity  arising  from  strong  pairing  in 
the  vicinity  of  the  Feshbach  resonance  [3,4].  Transition 
temperatures  of  up  to  half  the  Fermi  temperature  are  pre¬ 
dicted  to  result  from  the  strong  coupling  of  the  two-state 
Fermi  gas  to  the  bosonic  molecular  state  which  causes  the 
resonance.  Since  most  high  temperature  superconductors 
achieve  transition  temperatures  of  only  a  few  percent  of 
the  Fermi  temperature,  two-state  Fermi  gases  may  be  the 
highest  temperature  Fermi  superfluids  ever  studied  [5]. 
Further,  these  systems  may  permit  observation  of  the  tran¬ 
sition  from  weak  Bardeen-Cooper-Schrieffer  superfluidity 
to  Bose  condensation  of  strongly  bound  pairs  [6]. 

In  contrast  to  Bose-Einstein  condensates,  which  can  be 
prepared  and  studied  in  magnetic  traps,  two-component 
Fermi  superfluids  must  be  prepared  in  state-independent 
optical  dipole  traps,  since  the  required  pairs  of  hyperfine 
states  in  6Li  and  are  high-field  seeking  [1,2,7].  A  de¬ 
generate  Fermi  gas  has  been  produced  by  direct  evapo¬ 
ration  of  a  two-state  mixture  of  in  a  magnetic  trap, 
using  a  dual  radio-frequency -knife  method  [8],  Sympa¬ 
thetic  cooling  of  fermionic  6Li  to  degeneracy  also  has  been 
achieved  by  using  mixtures  of  6Li  with  bosonic  7Li  in  a 
magnetic  trap  [9,10].  However,  to  explore  superfluidity  in 
these  systems,  transfer  to  an  optical  trap  and  subsequent 
state  preparation  is  required.  The  procedure  for  prepar¬ 
ing  an  optically  trapped  two-state  Fermi  gas  can  be  greatly 
simplified  by  direct  evaporation  in  an  optical  trap. 

In  this  Letter,  we  demonstrate  all-optical  production  of 
a  degenerate  mixture  of  the  two  lowest  hyperfine  states  of 
fermionic  6Li  in  a  stable,  CO2  laser  trap  [11].  The  trap 
is  loaded  from  a  magneto-optical  trap  (MOT)  at  an  initial 


PACS  numbers:  05.30.Fk,  32.80.Pj 

temperature  of  150  jaK.  Degeneracy  is  obtained  by  forced 
evaporation,  accomplished  by  continuously  lowering  the 
trap  depth;  the  trap  is  then  adiabatically  recompressed  to 
full  depth.  At  this  stage,  more  than  105  atoms  remain 
at  temperatures  below  4  fi K,  less  than  half  of  the  Fermi 
temperature  of  8  /x,K.  TTiese  results  are  consistent  with 
scaling  laws  we  have  derived  for  the  phase-space  density 
as  a  function  of  trap  depth  [12]. 

Our  6Li  experiments  employ  a  C02  laser  trap  with  a 
single  focused  beam,  rather  than  a  crossed-beam  geome¬ 
try  as  used  recently  to  produce  a  87Rb  Bose-Einstein  con- 
denstate  (BEC)  by  forced  evaporation  [13].  Nevertheless, 
after  free  evaporation  at  full  trap  depth,  we  achieve  a  very 
high  initial  phase-space  density  of  —8  X  1CT3,  somewhat 
larger  than  that  obtained  after  free  evaporation  in  the 87  Rb 
BEC  experiments. 

A  commercial,  radio-frequency-excited  C02  laser 
(Coherent-DEOS  LC100-NV)  provides  140  W  at  A  - 
10.6  fim  for  the  trapping  laser  beam.  An  Agilent  (6573 A) 
power  supply  produces  stable  current  for  the  radio¬ 
frequency  source,  yielding  a  very  stable  laser  intensity. 
The  laser  output  is  deflected  by  an  acousto-optic  (A/O) 
modulator  to  control  the  power.  A  cylindrical  ZnSe 
telescope  corrects  the  output  of  the  A/O  for  ellipticity, 
and  the  beam  is  expanded  by  a  factor  of  10  before  passing 
through  an  aspherical  19.5  cm  focal  length  lens.  This 
lens  focuses  the  beam  into  the  vacuum  system,  yielding 
a  1/e2  intensity  radius  of  47  fi m.  The  corresponding 
Rayleigh  length  is  zq  —  660  fim.  With  an  incident  power 
of  65  W  in  the  trap  region,  the  trap  depth  is  estimated  to 
be  690  jjlK.  The  corresponding  radial  (axial)  oscillation 
frequency  for  6Li  is  predicted  to  be  6.6  kHz  (340  Hz), 
with  a  geometric  mean  of  v  =  (vxVyVt)x^  =  2400  Hz. 

The  radial  oscillation  frequency  is  measured  by  modu¬ 
lating  the  frequency  of  the  A/O  to  produce  a  sinusoidal  dis¬ 
placement  at  the  trap  focus  with  an  amplitude  of  0.2  /^m. 
After  the  sample  is  initially  prepared  at  a  temperature  of 
=*15  p, K,  the  modulation  is  applied  for  1  s.  The  number 
of  remaining  atoms  is  measured  by  resonance  fluorescence. 
Repeating  this  procedure  as  a  function  of  modulation  fre¬ 
quency  reveals  a  resonance  in  the  trap  loss  at  6.5  kHz,  in 
close  agreement  with  predictions.  Parametric  resonance 
methods  [14]  yield  results  consistent  with  the  expected 
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radial  and  axial  oscillation  frequencies  after  correction  for 
the  expected  resonance  frequency  shift  [15]. 

Extremely  low  residual  heating  rates  are  attained  in  the 
experiments.  At  the  maximum  trap  intensity  of  1.9  MW/ 
cm2,  the  optical  scattering  rate  is  2  photons  per  hour  as  a 
consequence  of  the  10.6  /im  wavelength  [16],  yielding  a 
recoil  heating  rate  of  only  16  pK/sec.  At  the  background 
pressure  of  <10“ 11  Torr,  heating  arising  from  diffractive 
background  gas  collisions  [17,18]  is  <5  nK/sec.  For  the 
trap  radial  oscillation  frequency  of  6.6  kHz,  the  intensity 
noise  heating  time  constant  is  estimated  to  be  >2.3  X 
104  sec  based  on  the  measured  laser  intensity  noise  power 
spectrum  [1 1].  A  residual  heating  rate  <5  nK/sec  is  mea¬ 
sured  at  full  trap  depth  over  200  sec.  Trap  \/e  lifetimes  of 
400  sec  are  observed. 

The  CO2  laser  trap  is  continuously  loaded  from  a  6Li 
MOT.  The  MOT  is  loaded  from  a  Zeeman  slower  for  5  sec, 
after  which  the  MOT  laser  beams  are  tuned  —6  MHz 
below  resonance  and  lowered  in  intensity  to  0.1  /sat  — 
0.25  mW/cm2  to  obtain  a  Doppler-limited  temperature  of 
—150  /iK  at  a  density  of  10ll/cm3.  Following  this  load¬ 
ing  stage,  the  MOT  gradient  magnets  are  extinguished  and 
the  upper  F  =  3/2  hyperfine  state  is  emptied  to  produce 
a  50-50  mixture  of  atoms  in  the  lower  |F  =  1/2, M  = 
±1/2)  states  [7]. 

The  |F  =  1/2,  M  =  ±1/2)  mixture  is  of  particular  in¬ 
terest,  as  it  is  predicted  to  exhibit  a  Feshbach  resonance 
near  850  G  [19].  A  convenient  feature  of  this  mixture  is 
that  the  .y-wave  scattering  length  vanishes  in  the  absence  of 
a  bias  magnetic  field  [19].  However,  the  scattering  length 
varies  between  0  and  — 300ao  as  the  bias  magnetic  field 
is  tuned  between  0  and  300  G  [19].  Hence,  rapid  evapo¬ 
ration  can  be  turned  on  and  off  simply  by  applying  or  not 
applying  a  bias  magnetic  field. 

The  number  of  trapped  atoms  is  enhanced  by  increasing 
the  intensity  of  the  CO2  laser  during  the  loading  stage 
[20].  To  accomplish  this,  the  beam  which  emerges  from 
the  trap  is  recollimated  after  a  ZnSe  exit  window  by  a 

19.5  cm  focal  length  ZnSe  lens,  and  then  retroreflected  and 
orthogonally  polarized  using  a  rooftop  mirror  oriented  at 
45°  to  the  incoming  polarization.  The  resulting  backward- 
propagating  beam  is  refocused  into  the  trap  region  through 
the  exit  lens.  After  passing  through  the  trap  region,  this 
beam  is  diverted  by  a  thin  film  polarizer  to  a  beam  dump  to 
avoid  feedback  into  the  laser.  Typically  1.5  X  106  atoms 
are  confined  in  the  forward  propagating  trap  beam  alone. 
The  backward-propagating  beam  increases  this  number  to 

3.5  X  106. 

After  the  CO2  laser  trap  is  loaded,  the  atoms  are 
precooled  by  free  evaporation.  To  initiate  evaporative 
cooling,  we  apply  a  bias  magnetic  field  of  130  G  by  revers¬ 
ing  the  current  in  one  of  the  MOT  gradient  coils,  yielding 
a  scattering  length  of  —  -100ao.  During  free  evapora¬ 
tion,  a  pneumatically  controlled  mirror  slowly  blocks  the 
backward-propagating  beam  by  diverting  the  power  into 
a  100  W  power  meter.  Since  this  beam  is  refocused,  the 
trap  region  is  Fourier-transform  related  to  the  plane  of 
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the  blocking  mirror,  and  the  trap  smoothly  evolves  into  a 
single  beam  configuration.  After  6  sec  of  free  evaporation, 
the  single  beam  trap  contains  N  —  1.3  X  106  atoms  at 
a  temperature  T  =  50  /iK.  This  precooling  procedure 
provides  excellent  initial  conditions  for  the  forced  evapo¬ 
ration  experiments,  since  the  resulting  phase-space  density 
for  each  state  at  full  trap  depth,  pi  =  (N/2)(hv)3/ 
(/cflT)3,  is  8  X  10“3,  which  is  extremely  high. 

In  all  of  our  experiments,  we  characterize  the  veloc¬ 
ity  distribution  of  the  trapped  gas  by  time-of-flight  imag¬ 
ing.  We  use  the  A/O  modulator  to  turn  off  the  CO2  laser 
trap  abruptly  (Ar  <  1  /is),  permitting  the  gas  to  expand 
for  a  time  between  400  /is  and  1.2  ms  in  zero  bias  mag¬ 
netic  field.  Residual  A/O  leakage  is  reduced  to  less  than 
2  X  10“5  of  the  maximum  intensity  by  extinguishing  the 
radio-frequency  synthesizer  output  prior  to  the  amplifier. 
Then  a  linearly  polarized  probe  laser  pulse  with  a  reso¬ 
nant  intensity  of  0.1/sat  and  a  detuning  of  3  half  linewidths 
(—9  MHz)  illuminates  the  gas  for  10  /is.  Simultaneously, 
a  noncopropagating  repumper,  resonant  with  the  D2  lines 
starting  from  the  F  =  3/2  state,  suppresses  optical  pump¬ 
ing  into  the  upper  F  =  3/2  hyperfine  state.  The  probe 
detuning  reduces  sensitivity  to  the  unresolved  excited  state 
hyperfine  structure  and  light  shifts  from  the  resonant  re- 
pumper.  For  the  selected  9  MHz  detuning,  the  expansion 
time  is  chosen  so  that  the  imaged  cloud  has  a  small  opti¬ 
cal  absorption  <35%.  An  achromat  at  the  vacuum  system 
exit  window  produces  a  1:1  image  of  the  atomic  distribu¬ 
tion  in  an  intermediate  plane.  This  plane  is  imaged  onto 
a  CCD  camera  (Andor)  using  a  microscope  objective  to 
produce  a  net  magnification  of  —4.  The  magnification  is 
calibrated  by  moving  the  axial  position  of  the  trap  focus 
through  ±1.25  mm  using  a  micrometer-controlled  transla¬ 
tion  stage.  Fitting  the  central  peak  of  the  distribution  to  a 
straight  line  yields  a  magnification  of  3.9. 

The  images  are  processed  to  obtain  the  transverse  spa¬ 
tial  distribution  by  integrating  the  measured  optical  depth 
in  the  axial  direction.  In  typical  measurements,  the  cloud 
expands  ballistically  by  100-200  /im  in  400  /is,  much 
larger  than  its  initial  transverse  dimension.  In  the  classical 
regime,  we  assume  ballistic  expansion  with  a  Maxwellian 
distribution.  In  this  case,  the  temperature  is  readily  deter¬ 
mined  from  the  transverse  \/e  width  of  the  cloud:  a(t)  = 
VyJ\/{2Trvr)2  +  f2,  where  v  =  j2kB7 '/M  is  the  ther¬ 
mal  velocity  and  t  is  the  time.  Since  vr  —  6.6  kHz,  for 
r  »  24  /is,  a{t)  =  vt.  Measurements  of  the  cloud  radius 
for  several  expansion  times  between  100  and  600  /is  fit 
very  well  to  a  straight  line. 

The  number  of  atoms  is  determined  from  the  spatially 
integrated  optical  depth  of  the  absorption  image  and  the 
absorption  cross  section  cr.  For  each  of  the  M  —  ±1/2 
magnetic  sublevels  of  the  populated  F  =  1/2  state,  a 
is  taken  to  be  (A2/7r)/[l  +  (2A/y)2],  2/3  of  that  of  the 
cycling  transition.  Since  the  excited  hyperfine  states  are 
unresolved  compared  to  the  linewidth  y  =  5.9  MHz,  this 
cross  section  contains  contributions  from  both  allowed 
transitions.  Results  for  the  number  are  consistent  within 
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10%  for  several  detunings  A  between  9  to  30  MHz 
and  -30  to  -9  MHz,  and  for  variation  of  the  camera 
focal  plane  over  ±  1  mm  from  the  plane  which  gives  the 
sharpest  image. 

After  precooling  by  free  evaporation,  further  cooling  is 
accomplished  by  lowering  the  trap  depth,  producing  forced 
evaporation.  We  have  developed  scaling  laws  for  the  num¬ 
ber  of  atoms,  collision  rate,  and  phase-space  density  as  a 
function  of  trap  depth  U  for  an  optical  trap  which  is  con¬ 
tinuously  lowered  [12].  These  scaling  laws  are  valid  for  a 
fixed  rj  =  U/{kBT)  »  1.  For  17  =  10,  we  find  that  the 
ratio  of  the  final  to  initial  phase-space  density  increases 
according  to  p/pi  =  (Ut/U)13.  This  result  shows  that 
lowering  the  trap  depth  by  a  factor  of  100  should  increase 
the  phase-space  density  by  a  factor  of  400,  producing  a 
degenerate  sample  for  pi  >  2.5  X  10“3.  To  maintain  a 
constant  value  of  77,  the  trap  should  be  lowered  from  its 
initial  depth  Ui  according  to  the  formula 

U(t)  =  Ui/{  1  +  t/rf,  (1) 

which  assures  us  that  the  lowering  rate  slows  as  the 
collision  rate  decreases  [12].  Taking  77  =  10,  we  have 
ft  =  1.45  and  1/r  —  2.0  X  10"3y,\  where  y,  is  the  ini¬ 
tial  elastic  collision  rate.  For  a  50-50  mixture  of  fermions, 
y.  =  ttNiM (Jv\ / (kBT()  with  Ni  the  initial  total  number 
of  atoms.  Note  that  y,  is  reduced  by  a  net  factor  of  4 
compared  to  a  single-component  Bose  gas  with  the  same 
parameters.  For  a  scattering  length  of  a  =*  -lOOao,  the 
elastic  cross  section  is  cr  =  87 to2  —  0.7  X  10“ 11  cm2. 
Using  Vi  =  (vxvyvz)l/3  =  2.4  kHz,  N(  =  1.0  X  106, 
and  Ti  —  50  pK,  we  obtain  y;  =  4.4  X  102  s”1  and 
r  =  1.1  sec. 

Unfortunately,  the  A/O  modulator  that  controls  the  CO2 
laser  intensity  produces  an  ellipticity  which  varies  as  the 
radio-frequency  (rf )  power  is  varied.  The  ellipticity  is  cor¬ 
rected  at  maximum  rf  power  by  a  cylindrical  telescope. 
However,  the  telescope  provides  only  fixed  compensation. 
Hence,  as  the  rf  power  is  decreased  to  lower  the  trap  depth, 
the  beam  becomes  elliptical,  reducing  vxvyvz  by  a  factor 
of  2  compared  to  that  expected  on  the  basis  of  the  laser 
power  alone.  Further,  we  find  that  the  direction  of  the 
beam  changes  by  3  mrad  as  the  rf  power  is  reduced  by  a 
factor  of  100,  causing  vignetting.  We  align  the  trap  beam  to 
minimize  this  vignetting,  but  beam  distortion  still  occurs. 
For  this  reason,  we  cannot  accurately  compare  our  evapo¬ 
ration  results  to  the  scaling  law  model.  To  compensate  for 
the  loss  of  confinement  arising  from  the  beam  distortion  as 
the  trap  is  lowered,  we  increase  r  to  3  sec.  The  trap  laser 
intensity  is  lowered  using  an  Agilent  (33120A)  arbitrary 
waveform  generator,  the  output  of  which  is  filtered  with  a 
time  constant  of  0.2  sec  before  being  applied  to  the  multi¬ 
plier  input  of  the  A/O  radio-frequency  generator. 

We  have  measured  atomic  velocity  distributions  after 
forced  evaporative  cooling  for  a  variable  time  tf.  To 
provide  a  calibrated  reference  trap,  time-of-flight  images 
are  recorded  after  adiabatic  recompression  to  full  trap 
depth  over  11  sec.  This  also  increases  the  spatial  den- 


FIG.  1 .  Absorption  image  (inset)  and  velocity  distribution  after 
10  sec  of  forced  evaporative  cooling  followed  by  recompression 
to  full  trap  depth.  An  average  of  five  trials  is  shown.  T /7>  — 
1,  as  determined  by  a  fit  to  a  Maxwellian  distribution.  vF  — 
20  cm/s  is  the  Fermi  velocity  for  a  total  N  —  8  X  103  atoms. 

sity  and  hence  the  Fermi  temperature,  while  preserving  the 
phase-space  density.  Figure  1  shows  the  velocity  distribu¬ 
tion  for  tf  =  10  sec.  The  total  number  of  atoms  remaining 
is  N  —  8  X  105,  corresponding  to  a  Fermi  temperature  of 
Tf  =  hv(6N/2)l^3/kB  —  15  p, K  for  each  state.  Assum¬ 
ing  a  Maxwellian  distribution,  the  gas  is  at  a  temperature 
of  15  /iK,  yielding  T/Tf  =  1.  At  this  temperature,  a  sub¬ 
stantial  number  of  atoms  have  velocities  greater  than  the 
Fermi  velocity  of  20  cm/sec. 

Near  degeneracy,  the  energy  of  the  atoms  contains  a  con¬ 
tribution  from  the  Fermi  energy  so  that  the  true  temperature 
is  lower  than  that  obtained  using  a  Maxwell-Boltzmann 
(MB)  distribution  which  assumes  that  all  of  the  energy  is 
thermal.  Hence,  the  low  temperature  absorption  images 
are  fit  using  a  Thomas-Fermi  (TF)  approximation  to  de¬ 
termine  T/Tf  [21,22],  where  the  Fermi  temperature  Tf  is 
calculated  using  the  measured  trap  frequencies  and  inte¬ 
grated  atom  number.  At  the  lowest  temperatures  achieved 
in  the  experiments,  the  MB  temperature  is  —10%  higher 
than  the  TF  approximation. 

Degeneracy  is  attained  for  tf  =  40  sec,  where  T  — 
5.8  p,K  and  T/TF  =  0.55  with  3  X  105  atoms  remain¬ 
ing.  At  this  temperature,  the  gas  is  degenerate,  and  p  “ 
(7>/T)3/6  —  1  [21].  We  have  also  measured  the  tempera¬ 
ture  of  the  atoms  in  the  lowered  trap  without  recompres¬ 
sion  to  full  trap  depth.  We  obtain  temperatures  a  factor  of 
—  10  lower,  i.e.,  —580  nK,  as  expected  for  a  harmonic  trap 
which  is  lower  in  depth  by  a  factor  of  —100. 

Figure  2  shows  the  velocity  distribution  for  //  — 
60  sec.  The  total  number  of  atoms  is  reduced  to  105, 
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FIG.  2.  Absorption  image  (inset)  and  velocity  distribution  after 
60  sec  of  forced  evaporative  cooling  followed  by  recompression 
to  full  trap  depth.  An  average  of  five  trials  is  shown.  T/TF  < 
0.5  as  determined  by  a  fit  using  a  Thomas-Fermi  approximation. 
vF  =  14  cm/s  is  the  Fermi  velocity  for  a  total  N  —  105  atoms. 


corresponding  to  a  Fermi  temperature  of  8  /xK.  The  mea¬ 
sured  temperature  is  below  4  /xK,  yielding  T/7>  =  0.48. 
Nearly  all  atoms  have  velocities  less  than  the  Fermi 
velocity  of  14  cm/sec. 

In  the  experiments,  we  achieve  high  evaporation  ef¬ 
ficiency  x  =  ^(pf/pi)/\n(Ni/Nf)  [23],  For  example, 
after  precooling,  but  prior  to  forced  evaporation,  Nj  ~ 
1.3  X  106  and  /?,-  =  8  X  10~3  per  state.  After  40  sec  of 
forced  evaporation,  Nf  =  0.3  X  106  and  p/  —  1.  Hence, 
X  ~  3.3.  The  overall  evaporation  efficiency  is  similar. 
Starting  with  the  loading  conditions  where  the  total  number 
of  atoms  is  N(  =  3,5  X  106  at  a  temperature  of  150  juK, 
we  obtain  x  ”  2.9  after  40  sec  of  forced  evaporation.  De¬ 
spite  the  trap  distortion  described  above,  these  results  are 
comparable  to  the  best  achieved  in  magnetic  traps  [23]. 

In  conclusion,  we  have  produced  a  degenerate,  two- 
component  6Li  Fermi  gas  in  a  single  beam  all-optical  trap 
by  direct  evaporative  cooling.  By  using  a  stable  CO2  laser 
trap  at  a  background  pressure  of  <10~u  Torr,  efficient 
evaporation  over  time  scales  of  85  sec  is  achieved.  In 
future  experiments,  it  will  be  possible  to  attain  scattering 
lengths  of  —  -300ao  by  increasing  the  bias  magnetic  field 
to  300  G,  thereby  increasing  the  elastic  cross  section  at  low 
temperature  by  nearly  a  factor  of  10.  This  should  enable 
preparation  of  a  degenerate  sample  in  just  a  few  seconds, 
producing  substantially  lower  temperatures  by  reducing 
the  detrimental  effects  of  any  residual  heating.  We  are 
currently  preparing  for  a  systematic  study  of  the  Feshbach 
resonance  at  higher  magnetic  field,  and  hope  to  observe 
superfluid  pairing  in  a  two-state  Fermi  gas. 
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We  derive  scaling  laws  for  the  number  of  atoms,  collision  rate,  and  phase-space  density  as  a  function  of  trap 
depth  for  evaporative  cooling  in  an  adiabadically  lowered  optical  trap.  The  results  are  in  excellent  agreement 
with  a  Boltzmann  equation  model  and  show  that  very  large  increases  in  phase-space  density  can  be  obtained 
without  excessive  slowing  of  the  evaporation  rate.  Predictions  are  in  reasonable  agreement  with  a  recent 
experiment  that  achieves  Bose-Einstein  condensation  by  evaporation  in  an  optical  trap.  We  also  discuss  evapo¬ 
ration  of  fermionic  mixtures  and  explain  why  Pauli  blocking  does  not  strongly  inhibit  cooling. 
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Far-off-resonance  optical  dipole  traps  are  well  known  to 
provide  nearly  state-independent  confining  potentials  for 
neutral  atoms  [1].  This  enables  study  of  systems  that  cannot 
be  stored  in  magnetic  traps,  such  as  stable  states  of  cesium 
atoms  or  diamagnetic  atoms  and  molecules.  Shallow  optical 
traps  have  been  used  to  confine  multiple  spin-state  spinor 
Bose  condensates  [2].  Optical  traps  also  are  likely  to  play  an 
important  role  in  studies  of  cold  neutral  fermions,  where  spe¬ 
cific  pairs  of  states  are  required  to  achieve  superfluidity 
[3-5].  A  degenerate  Fermi  gas  has  been  produced  by  direct 
evaporation  of  a  two-state  mixture  of  40K  in  a  magnetic  trap, 
using  a  dual  radio-frequency-knife  method  [6].  Sympathetic 
cooling  of  fermionic  6U  to  degeneracy  also  has  been 
achieved  by  using  mixtures  of  6U  with  bosonic  7Li  in  a 
magnetic  trap  [7,8],  However,  for  both  6Li  and  40K,  explo¬ 
ration  of  superfluidity  will  require  an  optical  trap.  Generally, 
the  procedure  for  producing  an  optically  trapped  degenerate 
gas  has  been  to  obtain  a  degenerate  gas  by  optical  cooling 
and  evaporation  in  a  magnetic  trap,  followed  by  transfer  to 
an  optical  trap.  However,  direct  production  of  a  degenerate 
gas  in  an  optical  trap  would  greatly  simplify  many  experi¬ 
ments,  and  has  been  explored  for  many  years  without  suc¬ 
cess  [9,10]. 

Recently,  a  Bose-Einstein  condensate  (BEC)  has  been 
produced  by  Barrett  et  al ,  using  forced  evaporation  of  87Rb 
in  a  C02  laser  trap  that  is  continuously  lowered  [11].  In  those 
experiments,  an  extremely  high  initial  spatial  density  is  ob¬ 
tained,  producing  a  very  large  elastic-scattering  rate  as  well 
as  a  high  initial  phase-space  density.  Hence,  the  trap  depth 
can  be  lowered  rapidly,  producing  a  BEC  in  a  few  seconds. 

In  this  paper,  we  describe  a  scaling  law  model  of  forced 
evaporation  in  a  continuously  lowered  optical  trap.  The  scal¬ 
ing  law  follows  from  a  simple  energy  evolution  equation  for 
the  trapped  atoms,  which  includes  the  energy  loss  arising 
both  from  evaporation  and  from  adiabatic  lowering  of  the 
trap  potential.  Previous  derivations  of  scaling  laws  for 
evaporation  in  magnetic  traps  have  not  explicitly  included 
the  time-dependent  potential  [12].  This  is  appropriate  for 
evaporation  in  traps  with  a  constant  potential  where  the 
radio-frequency-knife  method  is  employed  to  lower  the  trap 
threshold.  Since  the  trap  strength  does  not  change,  the  ex¬ 
plicit  time-dependence  of  the  potential  has  little  effect,  as 
shown  by  modeling  evaporation  from  a  magnetic  trap  with  a 
time-dependent  evaporation  threshold  [13].  To  determine  the 


correct  scaling  laws  for  adiabatically  lowered  optical  traps, 
where  the  trap  strength  varies  in  time,  we  explicitly  include 
the  time-dependent  potential  and  derive  scaling  laws  for  the 
phase-space  density,  number,  and  elastic  collision  rate  as  a 
function  of  the  well  depth.  Results  are  obtained  for  a  fixed 
large  ratio  tj  of  trap  depth  U  to  thermal  energy  kT ,  to  lowest 
order  in  the  small  parameter  exp(-  7),  using  a  harmonic  os¬ 
cillator  approximation.  For  77=  C//Xr7’  =10,  which  is  typical 
for  evaporation  in  optical  traps,  the  scaling  laws  show  that  by 
lowering  the  well  depth  by  a  factor  of  100,  the  phase-space 
density  can  be  increased  by  a  factor  of  400.  The  correspond¬ 
ing  elastic  collision  rate  is  reduced  by  a  factor  of  24.  By 
including  the  effects  of  loss  arising  from  background  gas 
collisions,  we  obtain  reasonable  agreement  with  the  results 
of  Barrett  et  al  [11]. 

It  is  well  known  that  the  evaporation  rate  of  a  gas  from  an 
optical  trap  of  fixed  depth  stagnates  as  the  temperature  drops 
[14].  At  low  temperatures,  the  number  of  colliding  pairs  of 
atoms  with  enough  energy  for  one  to  leave  the  trap  is  deter¬ 
mined  by  the  tail  of  the  Boltzmann  distribution.  Hence,  the 
evaporation  rate  is  suppressed  by  a  factor  Qxp(-UlkT).  For 
UlkT>  10,  the  evaporation  slows  dramatically,  and  it  is  nec¬ 
essary  to  force  evaporation  by  adiabatically  lowering  the  trap 
depth  [15]. 

The  optical  trapping  potential  can  be  written  generally  as 

t/(x,/)=-£/(/)g(x),  (1) 

where  g(x)  describes  the  trap  shape  and  g(|x|— with 
g( 0)  =  1 .  We  assume  that  evaporation  is  carried  out  at  low 
temperatures  near  stagnation,  where  the  average  thermal  en¬ 
ergy  kT<£U. 

To  determine  how  the  number  of  trapped  atoms  N ,  colli¬ 
sion  rate  y,  and  phase-space  density  p  scale  as  the  trap  depth 
U  is  lowered,  we  first  estimate  the  rate  of  energy  loss  from 
the  trap,  neglecting  atom  loss  arising  from  background  gas 
collisions.  Taking  the  zero  of  energy  to  be  at  the  bottom  of 
the  trap,  evaporating  atoms  will  have  an  average  energy  U 
+  akT ,  where  0=Sa<l  [14].  From  the  s-wave  Boltzmann 
equation  with  kT<U>  we  find  a=(7}-5)/(7}-4)  for  any 
potential  that  is  harmonic  near  the  minimum  [16].  The 
energy-loss  rate  arising  from  evaporation  is  then  N(U 
+  ctkT ),  where  N  is  the  rate  at  which  atoms  evaporate  from 
the  trap.  In  addition,  as  the  trap  depth  is  lowered  adiabati- 
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cally  at  a  rate  0>  an  energy  change  arises  from  the  change  in 
potential  energy.  Since  kT<U>  the  atoms  vibrate  near  the 
trap  bottom  in  an  approximately  harmonic  potential,  where 
EI2  is  the  average  potential  energy.  The  potential  energy  then 
changes  at  a  rate  (UIU)EI2  and  the  total  energy  E  obeys  the 
approximate  evolution  equation 

M{V+akT)+~.  (2) 

In  the  classical  limit,  E-ZNkT  is  the  total  energy  of  the 
trapped  gas  so  that  E=3NkT+3NkT.  Then,  the  contribu¬ 
tion  to  t  from  evaporation  is  proportional  to  N(U+akT 
-3 kT).  Hence,  the  cooling  rate  is  proportional  to  the  differ¬ 
ence  between  the  average  energy  carried  away  per  particle 
(=* U+akT)  and  the  average  thermal  energy  3 kT,  as  it 
should  be. 

Solving  Eq.  (2)  with  a  fixed  value  of  UfkT=  77,  the  num¬ 
ber  of  trapped  atoms  is  found  to  vary  with  trap  depth  as 

atU)  *  (3) 

where  i  denotes  the  initial  condition  at  /=0,  N=N(t)9  and 
U=U(t).  Here,  77'  =  77+0?=  77+  (77- 5 )/( 37 — 4). The  corre¬ 
sponding  phase-space  density  in  the  classical  regime  is  p 
=  N(h  v)*f(kT)39  where  v-  v(t)<xyfU  is  the  geometric  mean 
of  the  trap  oscillation  frequencies.  Using  Eq.  (3),  it  is  easy  to 
show  that  p=p(/)  scales  with  trap  depth  and  number  as 

Pi  UJ  UJ  • 


oms  such  as  6Li,  where  the  scattering  length  is  anomalously 
large,  evaporation  is  still  rapid  and  the  background  collision 
induced  loss  can  be  minimized  despite  the  reduced  collision 
rate. 

The  scaling  laws  are  derived  neglecting  background  gas 
collisions.  To  include  background  gas  collisions,  a  loss  rate 
term  -  T  bgE  is  included  in  Eq.  (2)  and  N  is  replaced  by 
Nevap>  where  the  evaporation  rate  Nevap~N+Y  bgN.  Then, 
one  can  show  that  the  scaling  laws  for  the  number,  collision 
rate,  and  phase-space  density  versus  trap  depth  are  reduced 
by  a  factor  exp(-T bgt)9  where  t  is  the  time  over  which  the 
trap  depth  is  lowered.  Different  exponential  factors  can  be 
derived  for  the  scaling  laws  versus  number.  Equation  (2)  is 
also  readily  modified  to  include  effects  of  residual  heating. 

Although  the  derivation  of  the  scaling  laws  does  not  ex¬ 
plicitly  include  the  time  dependence  of  the  trap  depth,  main¬ 
taining  a  constant  value  of  77  specifies  the  time  dependence 
£/(/),  which  follows  from  the  evaporation  rate.  To  lowest 
order  in  exp(—  77),  and  neglecting  background  gas  collisions, 
we  obtain  from  the  s-wave  Boltzmann  equation  [14], 

-2(77-4)  exp(-  77)7 N.  (6) 

Differentiating  Eq.  (3)  and  writing  y/y{  in  terms  of  UIUi9 
we  obtain 

m  1 

•  <7) 

where  the  time  constant  r  is  given  by 

12 

-=-v'(v-4)  exp( -  77) 7/ •  (8) 

T  J 


Equation  (4)  shows  that  for  77=  10,  p/pi=(C///C/)1,3.  For  an 
initial  phase-space  density  of  p,  =  3X10~3,  lowering  the 
well  depth  by  a  factor  of  84  yields  p=l.  For  an  energy- 
independent  scattering  cross  section,  the  elastic  collision  rate 
y=y(t)<*Nv*l(kT)  scales  with  trap  depth  as 


L 

7i 


*'/[2(*'-3)] 


(5) 


and  is  reduced  by  a  factor  of  21  for  a  factor  of  84  reduction 
in  well  depth  when  77=  10. 

An  important  feature  of  Eq.  (4)  is  that  the  increase  in 
phase-space  density  with  decreasing  number  is  identical  to 
that  obtained  using  a  radio-frequency-knife  method  with  a 
trap  depth  to  thermal  energy  ratio  of  77.  This  is  a  conse¬ 
quence  of  the  adiabatic  energy  loss,  which  ensures  that  the 
phase-space  density  does  not  change  as  the  trap  depth  is 
lowered  when  7^=0.  Hence,  the  phase-space  change  arises 
only  from  evaporation  as  in  the  radio-frequency-knife 
method.  For  an  optical  trap  with  77=10,  p/pj=(Af//JV)6,8, 
and  a  modest  decrease  in  number  leads  to  a  great  increase  in 
phase-space  density.  Unlike  evaporation  from  a  fixed  well, 
however,  the  collision  rate  for  an  energy-independent  elastic 
cross  section  decreases  as  ( UfUt )069  for  77=10,  and  run¬ 
away  evaporation  is  not  achieved  [17].  Nevertheless,  for  at¬ 


Including  background  gas  collisions  changes  / — ►[  1 
-exp(-r6g/)]/r/)g  in  Eq.  (7).  Equation  (8)  shows  that  the 
lowering  rate  scales  with  the  initial  collision  rate  7/ .  Accord¬ 
ing  to  Eq.  (7),  the  rate  of  decrease  of  the  well  depth  de¬ 
creases  with  time  as  the  collision  rate,  and  hence,  the  evapo¬ 
ration  rate  declines.  The  initial  elastic  collision  rate  for  a 
single  state  Bose  gas  in  a  harmonic  potential  is  y*(sec”  *) 
~47rNiM<rv*/(kTi)i  where  Nj  is  the  total  number  of  atoms 
initially  in  the  trap,  T,  the  initial  temperature,  cr=87ra2,  and 
vt  is  the  initial  mean  trap  oscillation  frequency  in  Hz.  For  a 
two-state  50-50  mixture  of  fermions  with  the  same  mass, 
scattering  length,  total  number,  and  trap  frequencies  as  in  the 
Bose  case,  the  rate  is  reduced  by  a  net  factor  of  4. 

The  scaling  law  predictions  are  valid  for  both  bosons  and 
fermions  in  the  classical  regime,  where  the  effects  of  quan¬ 
tum  statistics  can  be  neglected.  We  have  compared  the  scal¬ 
ing  law  predictions  to  a  detailed  Boltzmann  equation  model 
for  evaporation  of  a  two-state  mixture  of  fermionic  6 Li  in  a 
single  gaussian  beam  optical  trap.  Assuming  sufficient  ergod- 
icity,  the  evolution  of  a  low-temperature  trapped  gas  in  a 
time-dependent  potential  is  described  by  the  5-wave  Boltz¬ 
mann  equation  [13,14],  which  we  write  in  the  form  [18], 
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FIG.  1.  Phase-space  density  (each  state)  vs  trap  depth.  S-wave 
Boltzmann  equation  model  (dots);  scaling  law  predictions  (solid 
line).  The  time  varies  from  0  to  36  sec. 

The  left-hand  side  of  Eq.  (9)  describes  the  adiabatic  evolu¬ 
tion  of  the  gas  in  the  time-dependent  potential.  Physically,  in 
a  time  A t>  the  occupation  number  /(e,f)  changes  adiabati- 
cally  according  to  /(£,/  +  A/)  =/(6- A/(t?///<?/),0*  Hamil¬ 
ton’s  equations  require  (dHldt)  =  {U(€yt))y  where  the  angu¬ 
lar  brackets  denote  the  ergodic  average  of  the  time  rate  of 
change  of  the  potential  [18].  The  right-hand  side  of  Eq.  (9)  is 
a  Boltzmann  collision  integral,  which  redistributes  the  level 
occupation  numbers  f(e,t)  including  Fermi  statistics  [18]. 

A  comparison  of  the  results  of  the  scaling  laws  and  the 
Boltzmann  equation  predictions  for  a  two-state  mixture  of 
6Li  fermions  is  shown  in  Figs.  1  and  2.  The  scattering  length 
is  taken  to  be  a  =  -  300  a0,  which  can  be  obtained  by  apply¬ 
ing  a  magnetic  field  of  300  G.  The  trap  parameters  are  as 
follows:  Nt =  4  X  105  (total),  v,  =  1300  Hz,  7^=30  //.K,  U t 
=  300  /jlK,  i.e.,  77=  10  and  Tbg-  0.003  s"1.  The  well  depth 
is  lowered  according  to  Eq.  (7)  with  77=10  using  a  time 
constant  r  of  1.2  seconds,  as  predicted  by  Eq.  (8).  With  these 
parameters,  the  Boltzmann  model  yields  UlkT  between  10 
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FIG.  2.  Number  and  collision  rate  vs  trap  depth.  S-wave  Boltz¬ 
mann  equation  model:  Number  (triangles);  collision  rate  (squares); 
scaling  law  predictions  (lines).  The  time  varies  from  0  to  36  sec. 
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and  9.6.  Using  7 7=  10  in  the  scaling  laws,  we  obtain  nearly 
perfect  fits  to  the  Boltzmann  model  up  to  a  phase-space  den¬ 
sity  of  1,  where  the  effects  of  Fermi  statistics  become  impor¬ 
tant.  For  exponential  lowering  of  the  trap  depth,  where  77  is 
only  approximately  constant,  the  scaling  law  fits  are  in  rea¬ 
sonable  agreement  with  the  Boltzmann  model,  although  the 
fits  are  not  quite  as  good  as  for  constant  77. 

For  completeness,  we  describe  briefly  why  Fermi  statis¬ 
tics  does  not  severely  suppress  the  efficiency  of  evaporative 
cooling.  When  T<TF  the  collision  rate  within  the  trap  is 
reduced  to  Toe ycl(TITF)2<ych  where  ycl  is  the  classical 
collision  rate  evaluated  at  the  Fermi  surface,  T  is  the  tem¬ 
perature,  and  Tf  is  the  Fermi  temperature  [19].  The  factor 
(TITF)2  is  a  consequence  of  Pauli  blocking,  which  forbids 
collisions  into  occupied  energy  states,  as  observed  recently 
[20].  However,  in  evaporation,  one  of  the  final  states  is  es¬ 
sentially  unoccupied,  since  it  is  outside  the  trap.  Hence,  the 
evaporation  rate  is  suppressed  by  only  TITF ,  i.e., 

T  I  U—kTF\ 
re«a/*rc/-^r;exp|  kJ,  J. 

The  exponential  factor  describes  the  high-energy  tail  of  the 
Fermi  distribution,  which  is  responsible  for  evaporation 
when  the  trap  depth  U>kT.  This  is  essentially  the  same 
factor  that  appears  in  the  evaporation  of  a  classical  gas.  Since 
the  heat  capacity  also  scales  as  TITF  [21],  the  efficiency  of 
evaporative  cooling  in  lowering  the  temperature  is  not  seri¬ 
ously  compromised  in  a  two-component  Fermi  gas,  although 
the  sensitivity  to  residual  heating  is  increased.  Further,  the 
collision  rate  within  the  trap  is  always  fast 
compared  to  the  evaporation  rate  when  U—kTF>kT ,  since 
T/TF>  exp [-(U-kTF)/kT\  so  that  rethermalization  is  faster 
than  evaporation.  This  picture  explains  why  Pauli  blocking 
does  not  appear  to  strongly  affect  the  rate  of  decrease  of 
TITf  in  a  recent  theoretical  model  of  evaporation  for  a  two- 
component  Fermi  gas  [22]. 

It  is  interesting  to  compare  the  scaling  law  predictions 
with  the  experimental  BEC  results  obtained  by  Barrett  et  ai 
[11].  We  take  Nt  to  be  the  stagnation  value  after  1  second 
[11],  about  1/3  of  the  maximum  number  loaded.  In  the  actual 
experiments,  the  gas  is  not  given  time  to  stagnate  before  the 
trap  lowering  begins,  but  one  expects  that  most  atoms  are 
lost  in  a  small  fraction  of  a  second  before  the  well  depth 
changes  appreciably,  since  the  evaporation  slows  exponen¬ 
tially  as  the  temperature  drops  to  less  than  1/10  of  the  well 
depth  in  the  first  second.  For  the  final  conditions,  we  use  the 
data  for  a  laser  power  of  i>=350  mW,  near  the  transition 
between  the  classical  and  degenerate  regimes,  where  the 
scaling  law  is  approximately  valid.  We  use  the  trap  lifetime 
of  6  sec  and  a  lowering  time  of  2.5  sec  to  obtain  a  back¬ 
ground  loss  factor  of  exp(-F/,g0=exP(~2.5/6)  =  0.66.  The 
well  is  assumed  to  be  lowered  by  a  factor  UflUi 
-  ( vflVi)2=  1/84,  based  on  the  measured  trap  oscillation  fre¬ 
quencies.  We  take  7j=U/kT=\0  and  give  results  for  the 
case  where  all  atoms  are  in  a  single  hyperfine  state.  Table  I 
shows  that  the  predictions  are  in  good  agreement  with  the 
experiments  for  these  reasonable  assumptions  about  the  trap 
parameters. 
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TABLE  I.  Comparison  of  scaling  law  predictions  with  the  BEC 
experiment  of  Ref.  [11]. 


Initial  conditions 

Final  conditions 

Predictions 

Vj  =  1500  Hz 
N/=6.7XI05 

iy=  164  Hz 
w^i.sxio5 

Wy-l^XlO5 

r,  =  38  yuK 

7>=375  nK 

7>=450  nK 

y^^XlO3  sec"'1 

>y=300  sec-1 

>y=  372  sec*"1 

p,=  1/200 

Pf~  1.4 

ii 

In  conclusion,  we  have  derived  simple  scaling  laws  for 
the  number,  collision  rate,  and  phase-space  density  as  a  func¬ 
tion  of  trap  depth  for  atoms  in  a  time-dependent,  adiabati- 
cally  lowered  optical  trap  operating  near  stagnation.  The  in¬ 
clusion  of  the  time-dependent  potential  in  the  Boltzmann 
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evolution  equation  and  in  the  scaling  law  model  is  essential, 
as  it  accounts  for  the  adiabatic  energy  loss  in  the  trap,  which 
is  needed  to  preserve  the  phase-space  density  in  the  absence 
of  evaporation.  Our  results  show  that  the  phase-space  density 
in  optical  traps  increases  rapidly  as  the  trap  depth  is  lowered, 
when  the  ratio  of  trap  depth  to  temperature  is  large.  The 
reduction  of  the  elastic  collision  rate  with  well  depth  is  miti¬ 
gated  by  the  large  initial  spatial  [1 1,23]  and  phase-space  den¬ 
sities  obtainable  with  optical  traps,  as  well  as  the  large  scat¬ 
tering  lengths  obtainable  in  some  systems.  Hence, 
evaporation  in  optical  traps  appears  quite  promising  as  a 
means  to  achieve  degeneracy  in  a  variety  of  atomic  and  mo¬ 
lecular  systems.  Finally,  we  have  presented  a  physical  picture 
to  explain  why  Pauli  blocking  does  not  strongly  inhibit  cool¬ 
ing  by  evaporation  in  fermionic  mixtures. 
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